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g ■ ABSTRACT 

o ■ 

^ I We carry out a maximum-likelihood kinematic analysis of a sample of 1170 blue 

Q ' horizontal branch (BHB) stars from the Sloan Digital Sky Survey presented in Sirko 

■ et al. (2003) (Paper I). Monte Carlo simulations and resampling show that the results 

■ are robust to distance and velocity errors at least as large as the estimated errors 
. from Paper I. The best-fit velocities of the Sun (circular) and halo (rotational) are 
, 245.9ibl3.5kms~"'^ and 23.8ib20.1 kms~^ but are strongly covariant, so that f© — fhaio = 

> ; 222.1±7.7kms~^ If one adopts standard values for the local standard of rest and solar 
I motion, then the halo scarcely rotates. The velocity ellipsoid inferred for our sample is 

■ much more isotropic [(a^, a0,a^) = (101.4 ± 2.8, 97.7 ± 16.4, 107.4 ± 16.6) kms-^] than 

■ that of halo stars in the solar neighborhood, in agreement with a recent study of the 

■ distant halo by Sommer-Larsen et al. (1997). The line-of-sight velocity distribution 
. of the entire sample, corrected for the Sun's motion, is accurately gaussian with a 

'"q_|. dispersion of 101.6 ± 3.0kms~^. 

o : 

|H , Subject headings: Galaxy: halo — Galaxy: kinematics and dynamics — stars: horizon- 

^ I tal branch 

> ' 

tj ■ 1. Introduction 

P.: 

Among all visible components of the Galaxy, the stellar halo is the oldest and occupies the 
largest volume. It is therefore of interest both as a record of the early history of the Galaxy and 
as a tracer of the Galactic potential. 

The kinematics of halo stars in the solar neighborhood are well established, as shown in Table 1. 

For present purposes, the solar neighborhood extends as far as a few kpc from the Sun, but to a 
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distance less than the Sun's own Galactocentric radius. Evidently, the systemic rotation Vhaio of the 
local halo is small or zero, and the velocity ellipsoid is distinctly elongated in the radial direction. 

The more distant halo has been less well studied. Proper motions are generally not measurable 
for distant stars, so that line-of-sight velocities must be analyzed under a priori assumptions about 
the symmetries of the velocity ellipsoid. The stars are also fainter, reducing the accuracy of the 
measured quantities. 

Much work on the distant halo has been done using globular clusters, which have the advantage 
of being distinctive and bright, but have important drawbacks: there are only ~ 10^ of them, so 
that statistical uncertainties are large, and special conditions of formation or selective destruction 
by Galactic tides may have biased the properties of the cluster system compared to those of the 
field. Frenk & White (1980) analyzed the kinematics of a sample of 66 globular clusters under the 
assumption that the systemic rotation fhaio is constant; they concluded that the residual velocities 
are consistent with being isotropic within broad uncertainties. For a subsample of 44 metal-poor 
("F") clusters, they found Whaio = 42 it 33 km s~^ and a root-mean-square (rms) residual line-of-sight 
velocity uios = 124kms~^. In a similar analysis of clusters with [Fe/H] < —0.8, Armandroff (1989) 
found Vhaio = 46 lb 29 km s~^ and cios = 116 it lOkms"^. The restriction to low metallicities follows 
the discovery by Zinn (1985) that the more metal-rich clusters have a distribution and kinematics 
intermediate between those of the stellar halo and thin disk. 

Ratnatunga & Freeman (1989) sampled the halo field using relatively faint red giants at high 
Galactic latitude. Their summary for the more metal-poor stars near the south Galactic pole 
concludes that cr^ is approximately constant at ~ 75 kms~^ out to 25 kpc, which would be consistent 
with the solar neighborhood (Table 1) if the velocity ellipsoid is constant in cylindrical coordinates; 
however, for their most distant subsample near the SGP {{D) = 17.1 kpc), they measured o"ios = 
42± 10 (10 stars). 

In a study similar to that reported here, Sommer-Larsen et al. (1997) analyzed a sample of 
679 blue horizontal branch stars drawn from multiple sources and spanning Galactocentric radii 
7 — 65 kpc. With line-of-sight velocities and photometric distances, they found that the ra- 
dial velocity dispersion (in spherical coordinates) falls rather sharply from its value in the solar 
neighborhood, ar = 140 it lOkms"^, to an asymptotic value 89 it 19kms~^ at r > 20 kpc. By 
requiring dynamical equilibrium in an assumed logarithmic Galactic potential, they concluded that 
the tangential velocity dispersion should rise from its local value to erg ~ cr^ ~ 137kms~^ at large 
r. They did not claim, however, to have detected this rise directly in their data. They further 
suggested that the shift from radial to tangential anisotropy with increasing Galactocentric radius 
might support theories in which the halo formed by heirarchical merging of smaller systems. 

Because of the differences in kinematics cited above, it appears that local high- velocity, low- 
metallicity stars may not be entirely representative of the more distant halo. To study the more 
distant halo, we have isolated a sample of 1170 blue horizontal branch (BHB) stars from the Sloan 
Digital Sky Survey (Sirko et al. 2003, hereafter Paper I). In that paper, we estimated contamination 
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from foreground stars and blue stragglers to be less than 10% for bright stars {g < 18) and about 
25% for faint stars {g > 18). An improved color-absolute magnitude relation was also derived, 
and these absolute magnitudes are used throughout this paper to derive distance moduli. This 
sample of BHB stars provides effective tracers of the Galactic halo, with no kinematic bias. The 
sample may have selection bias, and is probably far from complete, but these biases do not affect 
the kinematic analyses presented in this paper other than degrading the statistics and increasing 
the errors. Using maximum-likelihood methods, we fit for the mean rotation and principle axes 
of the halo velocity ellipsoid on the assumption that these are constant in spherical or cylindrical 
coordinates. We also fit for the orbital velocity of the Sun. 

The present paper is organized as follows. We first discuss an overview of the kinematic analysis 
in §2 and appendix A. The results for the kinematic parameters of the halo are presented in §3 
and summarized in Table 2. Monte Carlo simulations which validate the analysis are presented in 
§4. §5 summarizes our main conclusions and discusses the agenda for future work. 



1.1. Terminology used in this paper 

The following terminology and conventions are used throughout this work. A heliocentric right- 
handed Cartesian coordinate system (x, y, z) is oriented so that +x points towards the Galactic 
center (GC) and +z points towards the North Galactic Pole, such that +y is the direction of disk 
rotation. Thus, = {vQ^x,'VQ,y,'VQ,z) = vlsr + {U,V,W) where {U,V,W) is the solar motion 
with respect to the local standard of rest (LSR). The velocity of the LSR is taken to be vlsr = 
(0, 220, 0) kms~^ (Kerr & Lynden-Bell 1986), and the canonical solar motion vector is {U, V, W) = 
(10, 5, 7) kms~^ (Dehnen & Binney 1998, rounded to integer values). Galactocentric spherical 
(r, 6, (f)) and cylindrical (i?, 0, z) coordinate systems will also be used. We make a rudimentary 
model of halo bulk rotation parametrized by a single velocity fhaio- A positive fhaio indicates a flat 
halo rotation curve in the same sense as the disk rotation. Radial velocities to stars are measured 
as redshifts which are subsequently and immediately corrected for the Earth's motion around the 
Sun and the Earth's rotation about its own axis. This heliocentric radial velocity Vr contains a 
component due to the solar velocity vq.^ The "line-of-sight" velocity v\os will be defined as Vr plus 
the projection of vq onto the line-of-sight to the star, minus the projection of its expected rotation 
component — fhalo^ (if any). That is, v\os is the component of the purely random component of 
velocity in a model of the halo with a flat rotation curve, projected onto the line of sight. We 
find v\os conceptually useful, but one should not forget that it is contingent upon our assumptions 
concerning the rotation of the Sun and the halo, whereas Vr is directly observed (modulo errors). 
Throughout, the distance from the Sun to the GC is assumed to be 8kpc (Reid 1993). 



Here the subscript "r" is used for components along the Une of sight, but elsewhere it will be used for the radial 
component in spherical coordinates based on the Galactic Center. The context should make our meaning clear. 
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2. The velocity ellipsoid of the Galactic halo: overview 

The velocity ellipsoid is defined such that its axes align with the principal axes of the velocity 
dispersion tensor, and the lengths of its axes cr are equal to the components of the velocity dispersion 
along those directions. It is difficult to know the orientation of the velocity ellipsoid so it is 
reasonable to assume it is oriented along galactocentric spherical or cylindrical coordinates at every 
point in space (Binney & Merrifield 1998). Even with ~ 10^ stars, it is difficult to assess whether 
and how the components of the velocity ellipsoid vary in space, so they are assumed constant, 
although we test this below by dividing our sample into a few magnitude bins. 

A maximum likelihood technique is used, as outlined in Appendix A. A uniform prior is 
adopted (Lupton 1993): when the probability of observing the set of stars at the observed positions 
and heliocentric radial velocities given an assumed cr is maximized as a function of <t, the probability 
of <T being the true velocity ellipsoid is maximized. In symbols, 0(uios I "") oc (j){cr \vios) where (p 
represents probability density. Naturally, (j){v\os \ a) depends implicitly upon the orientation of the 
velocity ellipsoid with respect to the line of sight, and hence upon the angular position {a, 6) and 
distance of each star. It is straightfoward to expand the maximum likelihood procedure to include 
four additional parameters: the three components of the velocity of the Sun with respect to the 
galactic center, v©, and the bulk rotation velocity of the halo, fhaio- 

We assume that the bulk rotation of the halo corresponds to a uniform velocity in the direction 
of galactic rotation. It is difficult to apply this analysis on this sample to more complex halo rotation 
laws because of the moderate degeneracy between the halo rotation and VQ^y (see figure 2, discussed 
below); this degeneracy would increase if one attempted to fit for more parameters of the bulk halo 
rotation. 

Distance moduli are derived from the color-absolute magnitude relation of Paper I. The mea- 
sured angular coordinates of the stars {a, 6) are virtually error-free. However, the magnitudes, and 
thus distances, are susceptible to several sources of error: extinction, photon statistics, and intrin- 
sic brightness variations from star to star which create deviations in the stars' distance moduli. 
Therefore, magnitude errors should be included in the likelihood analysis. As shown in Paper I, 
typical magnitude errors are about Am 0.2. In anticipation of our finding, discussed below in 
§3, that the kinematic analysis is very robust with respect to magnitude errors, we adopt an even 
more conservative magnitude error: Am = 0.5, and assume the errors follow a gaussian distri- 
bution. Equation A3 shows the probability density <j) of observing a star at a given position in 
space, line of sight velocity vios, and cr. Having assumed that magnitude errors have a gaussian 
distribution with some width Am, the likelihood of observing a star considering magnitude errors 
can be found by integrating the probability density of equation A3, multiplied by a gaussian weight 
function for the magnitude error, along the line of sight. Figure 1 shows an illustration of this 
procedure for two stars with artificial data. These two stars are constructed to have (/, 6) = (0, 30), 
g = 15.5, and absolute magnitude 0.7, which places them approximately 9 kpc from the Sun above 
the galactic center. The evaluation of the probability density (j) for each star alone along the line 
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of sight is shown in the upper right panel. One star's hne-of-sight velocity is vi^g = Okms~^ (solid 
line and filled dots) and the other's is vi^s = lOOkms"^ (dotted line and open circles). We apply 
<T = (150,50, 100)kms~^ (spherical coordinates) to this model (note that the cp component of the 
velocity ellipsoid is irrelevant for this particular case). It is easy to see several trends in the upper 
right plot of figure 1. If there were no distance measurement, for example, the stationary star would 
be most likely to be found at ~ 7kpc because the line of sight coincides there with the smallest 
(9) axis of the velocity ellipsoid. In fact, the probability density of observing a stationary star at 
that point is just {27rag)~^^^ ^ .008skm~^, as shown. Also, the receding star is least likely to be 
observed there because its viqs can most easily be "explained" by most of its contribution coming 
from the radial component of the velocity ellipsoid. Also note that the stationary star always has 
a higher probability than the receding star; this follows from the gaussian nature of the probability 
of observing vios in equation A3. 

To incorporate the magnitude error Am into the analysis, the probability density as a 
function of distance must be multiplied by a gaussian centered on the measured magnitude with 
width Am, and then integrated over distance. This integration is most easily accomplished with 
Gauss- Hermite integration, which samples (j) at A'^p = 11 points, shown in the upper right panel 
of figure 1 as dots and circles. These 11 values of (j) at the specified abscissae are multiplied by 
the coresponding weight factors Wj (Abramowitz & Stegun 1972), which are plotted in a discrete 
fashion as the lower line in the middle panel of figure 1. The consideration of magnitude errors also 
introduces Malmquist bias. If the number density of observable BHB stars in the Galactic halo 
follows a spherical power law p oc r"° in Galactocentric radius, then Malmquist factor M(r') oc 
^/3j,-o j]2ust weight the probability density (j), in addition to the gaussian weight Wj due to the 
magnitude error, where r' is heliocentric distance. These weighting factors are normalized so that 
"^f^ WjM{r) = 1. The function M{r') is shown in the middle panel of figure 1 as the top line, 
assuming a = 3.5. Note that for a = 3, M[r ^ 8kpc) constant. 

Finally, the probability density of observing a star at a given angular position, distance mod- 
ulus, and line-of-sight velocity, with assumptions on cr. Am and a for the Malmquist term, is just 
(j){r)wjM{r). The term cj){r)wjM[r) is plotted in the lower panel of figure 1 for our two fictitious 
stars; again the stationary star is the solid line and the receding star is the dotted line. When these 
terms are summed over j, we arrive at new determinations for the probability density of observing 
the star, indicated in the upper right panel of figure 1 by the horizontal lines. The original estimate 
of (f) was given by the middle dot and circle (at ~ 9kpc), so it is evident that the "hill" and the 
"dip" in (p{r) for the two stars are responsible for modifying the original estimate of 4> up or down 
in the correct direction. 

The quantity a for the halo is usually taken to be between about 3 and 3.5 (Yanny & Newberg 
et al. 2000, and references therein). In the kinematic analysis below we assume a = 3.5, following 
the study of horizontal branch stars of Kinman, Suntzeff, &: Kraft (1994). However, we find our 
results are fairly insensitive to a; even a = 10 results only in small changes to our final results. 
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Another possible source of measurement error is in the radial velocities. As determined in 
Paper I, typical velocity errors may be as high as ~ 30kms~^. According to equation A3, the 
probability density of observing a star at a certain line-of-sight velocity v\os (with respect to the 
GC) follows a gaussian distribution, centered on vios = Okms~^, with a dispersion aios that depends 
on geometry and cr. Thus the probability density of observing the star at viq^ considering a possible 
velocity error that follows a gaussian distribution with dispersion Av is just 



note this reduces to the original equation A3 as Av — > 0. 

We have discussed the various modifications available to the maximum likelihood technique. 
In summary, they are the following: 

• Subsample. Although we first determine the velocity ellipsoid using all stars, we can extract 
subsamples, specifically in certain magnitude bins, to investigate whether the velocity ellipsoid 
is a function of position in space. 

• Spherical or cylindrical coordinates. The choice affects the values of x in the maximum 
likelihood procedure but is otherwise straightfoward. 

• Magnitude errors. We either ignore magnitude errors or else use 11-point Gauss-Hermite 
integration with Am = 0.5 and a = 3.5. 

• Velocity errors. We either ignore velocity errors or else use equation 1 with Av = 30kms~^. 

• Number of fitting parameters. Several combinations of fitting parameters are used. We fit 
the three cr parameters keeping (v0, fhaio) at canonical values; we fit (fhaio, cr) keeping at 
its canonical values; and we also fit all seven parameters simultaneously. 



Table 2 contains the determinations of (v©, Vhaioi ^) for several combinations of model assump- 
tions discussed above; the first four columns describe these model assumptions. The first column 
states what criteria were used to select subsets of BHB stars from the total set of 1170 BHB stars, 
and the number of stars in each subset is listed in the second column. The third column indicates 
whether a cylindrical or spherical geometry was assumed. The fourth column indicates how errors 
were handled; "0" indicates that errors were not considered. As described in §2, "m" indicates that 
11-point Gauss-Hermite integration was used with Am = 0.5, and "v" indicates that velocity errors 
were considered with At; = 30. The following columns show the parameter determinations. All 
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determinations in table 2 are for the three-parameter fit with (v0,Uhaio) fixed at canonical values 
(Dehnen & Binney 1998) as shown, except for the second and third lines which fit four and seven 
parameters, respectively. The quoted errors consider correlations among the parameters, and the 
errors in parentheses do not consider the correlations. The former error is defined as —Wi/2 and 
the latter as 1/ ^J —2Mi where Mj represents the diagonal elements of the second-derivative matrix 
of the log-likelihood function, and Wi represents the diagonal elements of the inverse of that matrix. 

The first row in Table 2 shows the results for the velocity ellipsoid for the complete sample 
(-^stars = 1170) in spherical coordinates, with no magnitude or velocity errors considered. The Sun's 
velocity was fixed at the canonical value of vq = vlsr + (C/, F, W) = (0, 220, 0) -|- (10, 5, 7) kms"-*^ 
(Dehnen &: Binney 1998). Under these assumptions all three components of the velocity ellipsoid are 
about lOOkms"^. The second row in Table 2 shows the results of simultaneously fitting fhaio and 
to the same sample and otherwise same set of assumptions discussed above. The velocity ellipsoid 
has essentially the same determined values, and the rotation velocity of the halo is consistent with 
zero. The third row in Table 2 is the seven parameter fit to the sample. Again, the velocity ellipsoid 
is essentially the same, and the rotation velocity of the halo is seen to be marginally consistent with 
zero. Reassuringly, the velocity of the Sun with respect to the galactic halo in the y direction (or 
equivalently, the Sun's velocity in an inertial frame assuming the halo is actually non-rotating) is 
close to the expected value: VQ^y — Uhaio = 222.1 ± 7.7kms~^. The error follows from rotating the 
eigenvectors by 45 deg and is smaller than either error for v^^y or f halo because the two quantities 
are strongly covariant. The x and z directions of Vq are new measurements of the solar motion 
{U,W). Our value of U is consistent with that of Dehnen &; Binney (1998); however our value of 
W differs by several a. 

Figure 2 shows confidence regions for the components of the velocity ellipsoid for this case. 
Each panel in figure 2 is a two-dimensional parameter space where the other five parameters are 
held at their determined values and assumed to have symmetric errors, so the confidence regions 
are only approximatations to those one would obtain by integrating over the other five dimensions. 
Note that (Jr is very well constrained. This is because the typical distance to the halo stars in our 
sample is large compared to the distance from the Sun to the galactic center, so measurements of 
v\os from Earth very nearly align with direct measurements of the galactocentric radial component 
of velocity. Also, each component of the velocity ellipsoid is anticorrelated with the other two; this 
demonstrates an intuitive "conservation law" of velocities. In the lower left panel, it is apparent 
that f halo and v^^y are correlated, so that the quantity VQ^y — f halo is more tightly constrained than 
either quantity alone. This of course follows from the reasonable assumption that the halo rotation 
at the location of the Sun is in the y direction. 

The result that ar ~ 100 km s^"*^ contrasts with most estimates of the velocity ellipsoid based 
on stellar samples in the solar neighborhood (WooUey 1978; Pier 1984; Carney & Latham 1986; 
Norris 1986; Layden et al. 1996). We now investigate some other sets of assumptions that may 
influence the results. 
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The fourth hne in Table 2 introduces magnitude errors into the set of assumptions. These 
errors, at least in our implementation of a gaussian with width Am = 0.5, apparently have little 
effect on the final results. Note that for a given kinematic model, distance affects the distribution 
of I'los only through the orientation of the velocity ellipsoid with respect to the line of sight. If the 
dispersions were isotropic and the halo nonrotating, which seem to be nearly the case, then vyos 
would be independent of distance. This may explain the insensitivity to magnitude errors. The 
fifth line in Table 2 incorporates both magnitude and velocity errors. Again, the final results are 
seen to remain relatively constant; however the values of the velocity ellipsoid components decrease 
modestly. In light of equation A3, the probability of observing a star at a certain velocity vy^s 
is always higher for smaller values of |uios|; thus velocity errors tend to "favor" smaller velocities, 
which propagate to the smaller values of the velocity ellipsoid. However, the effect is small, of order 
Vo"^ + Av'^ — (7 ~ Skms^"*^. 

The next three lines in Table 2 show three-parameter determinations for the complete sample 
of BHB stars assuming the velocity ellipsoid is oriented along cylindrical coordinates. As for the 
spherical coordinates case, all three components of the velocity ellipsoid are still about lOOkms"^. 

Determinations for cylindrical scenarios incorporating magnitude and velocity errors are shown 
in the tables, and as before, do not appear to affect the final results appreciably. Since velocity 
errors behave in a predictable way and have little effect on our final results, and magnitude errors 
have no noticable effect, we do not consider these errors further. It should be remembered that our 
error estimate was based on a set of highly simplified, and perhaps dubious, assumptions that do 
not take into account systematic errors, for instance. 

In the following rows in Table 2 the sample of stars is modified. There are 733 stars with g < 18, 
so these stars extend from ~ 7 kpc to ~ 30 kpc. The sample with < 16 contains the 227 brightest 
stars (distances out to 11.5 kpc). The velocity ellipsoid is seen to remain remarkably constant down 
to distances from the Sun only marginally outside the coverage of previous surveys (e.g., V ^ 13 
for the sample of Layden et al. (1996)). The most intriguing result is that the radial component 
of the velocity ellipsoid ar (or aji in cylindrical coordinates) just outside the solar neighborhood is 
significantly smaller than the canonical value of ~ 150kms~^ in the solar neighborhood (Table 1). 
This incompatibility will be further discussed in §5. 

The sample with g > 18 contains the 437 faintest stars. Since the stars are all very far away, the 
non-radial components of the velocity ellipsoid are very poorly constrained. Note that the quoted 
errors may also be ill-determined because of departures from gaussianity for small eigenvalues. 
However, the essential result remains the same: the radial component of the velocity ellipsoid is 
~ 100 km s-^ 

Finally, to address concerns that the determination of the velocity ellipsoid may be contam- 
inated by halo substructure, the last two lines of Table 2 show the results for all stars excluding 
the Sagittarius stream. This exclusion is conservatively accomplished by simply rejecting all BHB 
stars within 10 degrees of the great circle with pole {i, b) = (94, 11), which approximates the Sagit- 
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tarius stream (Ivezic et al. 2003, Paper I). Since the Sagittarius stream is the most conspicuous 
substructure in our sample of BHB stars (Paper I), the stabihty of the velocity ellipsoid solution 
indicates that halo substructure does not appreciably influence the results. 

4. Monte Carlo simulations 

An independent way of verifying the accuracy of the maximum likelihood technique is through 
Monte Carlo simulations. In case the geometry of the SDSS footprint (see Paper I) affects the 
analysis, we assemble 1170 artificial stars with the same angular coordinates and measured g 
magnitudes as the complete BHB sample, but with newly constructed velocities. For simplicity 
absolute magnitudes are assumed constant at 0.7. Only velocity ellipsoids aligned with spherical 
coordinates will be considered here. These Monte Carlo simulations perform the three-parameter 
fit to the velocity ellipsoid with = (10, 225, 7) kms~^ and I'haio = Okms~^. Note, however, that 
the values of vq and fhaio are irrelevant here since they are first added to the stars' velocities in 
construction, and then subtracted as known in the parameter determination. 

We run 100 Monte Carlo simulations on each of four setups. The first setup does not consider 
magnitude or velocity errors; the second setup allows for Am = 0.5; the third assumes Av = 
30kms~^. The fourth setup considers the magnitude and velocity errors as well as allowing for 
contamination. 

The constructed set of artificial stars has known space positions, so it is straightforward to 
apply an assumed velocity ellipsoid to obtain a simulated three-dimensional space velocity for each 
star. Here we assume <t = (150, 100, 80) km s-^ to strain the maximum likelihood analysis. The 
three-dimensional velocities are then projected onto the line of sight to each star to obtain viq^, and 
then the projections of vq and fhaio are subtracted to obtain Vr- This comprises the simulated data 
set for the first setup. To incorporate magnitude errors, the radial distances to the stars are shuffled, 
i.e. the g magnitudes are shifted by a gaussian error distribution of width Am, Vy is calculated as 
before, and then the g magnitudes are shifted back to their "observed" positions. Velocity errors 
are incorporated by adding an error to each observed Vr- To incorporate contamination, a fraction 
of stars are randomly flagged as disk stars. These stars are then assumed to be 2.5 times closer 
(2 mag) to the Sun than they would be if they were BHB stars. Their velocities are taken to 
be 220kms~^ in the — (/> direction plus three random components from a disk velocity dispersion, 
which is assumed isotropic with components (20, 20, 20) km s~^. Dehnen & Binney (1998) find a 
non- isotropic disk velocity dispersion with components of the same order (depending on color), but 
our intent is to contaminate the simulated halo star population in a simple and controlled way. 

Figure 3 shows the results of these Monte Carlo simulations with the four setups colored black 
(no errors), blue (magnitude errors), red (velocity errors), and green (magnitude and velocity errors 
with 10% contamination). As shown, the scatter in the results of the Monte Carlo simulations is 
comparable to the confidence regions shown in figure 2, so our error analysis is verified to be 
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reasonable. As in the analysis of the real data, the radial component of the velocity ellipsoid ar is 
very well constrained; one can also see anticorrelation between any two components. The black dots 
and blue dots occupy essentially the same region of parameter space; this verifies, independently of 
the method used in §3, that magnitude errors of the order Am do not affect our results significantly. 

As shown in Table 3, the mean determinations of a over 100 Monte Carlo simulations recover 
the input values for the "no errors" case and the "magnitude errors only" case. It is also clear that 
velocity errors of scale Av = 30kms~^ push the measured values higher by ~ 5kms~^ because the 
error adds in quadrature. This result is consistent with the independent determination in §3 which 
stated that with Av = 30kms~^ the "true" values of the velocity ellipsoid (i.e., those measured 
using the analysis that allows for velocity errors) were ~ 5kms~^ lower than the "measured" 
values (not considering velocity errors). Finally, it is interesting that for the case which includes 
contamination by disk stars, the determination of a^j) increases while those of ar and ae both 
decrease. This is because the disk stars are constructed to travel primarily in the (j) direction. 

However, the main conclusion to be gained from these Monte Carlo simulations is that the 
radial component of the velocity ellipsoid never strays too far from its correct value of 150kms~^, 
even allowing for errors such as contamination. Thus the maximum likelihood analysis is verified 
to be reliable and stable. 

5. Summary and discussion 

Applying maximum-likelihood methods to the line-of-sight velocities and photometric distances 
of BHB stars, we have modeled the kinematics of the relatively distant galactic halo. (The median 
distance of the sample is about 25kpc.) No significant anisotropy of the velocity ellipsoid is de- 
tected. In a model aligned with spherical coordinates, for example, a = (101.4,97.7, 107.4) kms~^. 
However, the errors in ag and are much larger, ~ 15 — 20kms~^, than that of ct^, because of 
the great distances of most stars: that is, the Galactocentric radius r of most stars is more than 
twice that of the Sun, so that the direction of the line of sight is not far from radial despite the 
broad angular range of our sample. This result contrasts starkly and surprisingly with the radial 
anisotropy of halo velocities in the solar neighborhood (cr ~ (150, 100, 100) km s~^; see Table 1). 
Our results are consistent, however, with the only other comparable study of distant halo field 
stars: Sommer-Larsen et al. (1997). 

We are able to fit separately for the rotation speed of the halo (fhaio) find the orbital velocity 
of the Sun, although the two are strongly covariant. The halo rotation is marginally consistent 
with zero at the one-sigma level, which is ~ 20kms~^. The difference VQ^y — t^haio = 222.1 it 
7.7kms~^ is more tightly constrained than either quantity separately. This is consistent with a 
non-rotating halo if one accepts the lAU recommendation for the local standard of rest, vlsr = 
220kms~^ (Kerr & Lynden-Bell 1986) and Dehnen & Binney (1998)'s solar motion, V0 — vlsr = 
(10, 5, 7) kms~^ , although some authors have advocated a substantially smaller (Oiling & Merrifield 
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1998) or larger (Uemura et al. 2000) value of the azimuthal component of vlsr- (If the solar rotation 
is substantially less than 225kms~^, then according to our results, the halo should counterrotate!) 
Our result is also statistically consistent with most of the studies cited in Table 1, although the 
latter seem to prefer a slight positive halo rotation. Note that what those studies actually measure 
is the mean velocity of local halo stars relative to the Sun; the values for Vhaio in the fourth column 
of that table have been adjusted to Vq = 225kms~^ so that they can be compared directly to 
our result. Presumably, vlsr will soon be determined with exquisite accuracy by a variant of the 
classical visual-binary method applied to stars orbiting Sgr A*, but accurate values are not in the 
public domain at the time of writing (Ghez et al. 2003; Salim & Gould 1999). 

Of course, the measurement of Uhaio depends upon the assumption that the halo rotation curve 
is flat; line-of-sight velocities would be unaffected by adding a solid-body component to the angular 
velocity of the halo, although proper motions would be (see below). 

These are the main results of our study. It is likely that more science could be extracted from 
this sample, and completion of the SDSS will allow it to be enlarged. For the present, we simply 
highlight certain interesting features of the data that have not yet been thoroughly explored. 

The distribution of vios is plotted in figure 4 and beautifully follows a gaussian of width 
cjios = 101.6 lb 3.0 kms^^. The uncertainty quoted here is one-sigma statistical, to which should 
be added a similar systematic uncertainty pending better quantification of our measurement er- 



rors. An estimate of this systematic error is easily determined as crios — \/ cr? — Av'^ ~ 3.4 km s 



{Av = 26kms~^ from Paper I). The line-of-sight velocities (figure 4) are about as gaussianly dis- 
tributed as the number of stars would allow. This is consistent with (although not required by) the 
isotropy of the velocity ellipsoid that we find, since a strongly anisotropic ellipsoid, even if gaussian 
along each individual axis, would not produce an exact gaussian in an average over many lines of 
sight. Gaussianity could result from large observational errors, regardless of the intrinsic velocity 
distribution, if the errors were gaussian, but our errors seem not to be large enough to dominate the 
observed distribution (Paper I). By contrast, the distribution of transverse velocities of halo stars 
in the solar neighborhood is distinctly nongaussian, with a negative kurtosis (Popowski & Gould 
1998). This is another indication that the local halo differs from the more distant halo probed by 
our BHB sample. 

Another aspect of the data that could be thoroughly explored is the smoothness of the BHB 
phase-space distribution. The notion of a stellar velocity ellipsoid implicitly presumes that phase 
space is smoothly populated. If instead the distribution is very lumpy and the number of important 
lumps or streams is small, then statistical characterizations of individual halos, in particular that of 
the Galaxy, may not be meaningful. In cosmology, however, statistical descriptions are unavoidable. 
Current theory holds that galactic halos have been built up by accretion of lesser systems (Searle 
& Zinn 1978; Bullock, Kravtsov, & Weinberg 2001). Accreted gas joins galactic disks, and the 
details of its history are effaced by dissipation. Stars that form in these systems before they are 
accreted (and presumably also their dark matter), however, behave collisionlessly and preserve a 
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partial record of their origins. The structure of the stehar and dark-matter halos should be related. 
Helmi, White, & Springel (2003) conclude from a detailed analysis and scaling of cosmological 
simulations that the number of streams composing a halo similar to that of the Galaxy should be 
large, amounting to almost 10^ within 25kpc, and that most streams joined the Galaxy early on, 
allowing more time for phase mixing. The number of our BHB stars is much smaller than this, so 
that one might expect to find a smooth distribution. And yet we see structure in our sample that 
is consistent with the tidal stream of the Sagittarius Dwarf (Paper I), as previously delineated in 
SDSS data (Ivezic et al. 2000; Yanny k Newberg et al. 2000; Ibata et al. 2001a; Newberg k Yanny 
et al. 2002) and in a survey of carbon stars (Ibata et al. 2001b). Most of our results for the halo 
velocity ellipsoid are scarcely affected by excluding that part of the sample that lies close to the 
Sagittarius stream on the sky. In short, it is not yet clear whether the kinematics of the Galactic 
halo admit a statistical decription. Larger samples of stars will be necessary to clarify this. 

In Paper I, we used proper motions only as a check on the purity of our sample. Calibration 
against the QSOs shows that the proper motions of individual BHB stars in our sample are smaller 
than measurement errors by factors ~ 3 — 5, but statistical detection of the kinematics of the 
population as a whole should certainly be possible. Among other benefits, this would allow a 
direct constraint on the solid-body component of halo rotation and a further test of the isotropy of 
the velocity dispersions. Systematic errors in the proper motion data present many complexities, 
however, so we defer their analysis to a later paper. The quality of the USNO-B data should be 
appreciably better than that of the USNO-A, and one should probably use the former to determine 
the proper motions of SDSS objects. 

Our BHB radial velocity data could be used to constrain the mass of the outer Galaxy. We have 
not yet attempted a serious dynamical analysis, partly in order to do full justice to the kinematics, 
but also because mass measurements rely more heavily on aspects of the data that are less certain: 
namely, the photometric distances and the completeness as a function of magnitude. However, since 
no kinematic information was used in the selection of the sample, we may assume that the velocities 
fairly represent the whole population of BHB stars, whose spatial density is often described as a 
powerlaw in Galactocentric radius, n tx r~", with exponent a ~ 3.5 (Kinman, Suntzeff, k Kraft 
1994). In that case, to the extent that the dispersions are isotropic and constant with radius and 
rotation is negligible, and that the potential is spherical, Jeans' equations imply 



The value of Vchc so obtained, with a ~ lOOkms" , is substantially lower than the circular ve- 
locity in the vicinity of the Sun. Faced with similar evidence drawn from their own independent 
BHB sample, Sommer-Larsen et al. (1997) concluded that the tangential velocity dispersion of 
the halo probably rises beyond r ~ 20kpc so as to exceed the radial dispersion and be consistent 
with a constant circular velocity. At such distances, however, the tangential dispersions are very 




(2) 



(3) 
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poorly constrained by line-of-sight velocities without a priori assumptions about the potential (see 
Table 2). Rather than make such assumptions, one may prefer to await further direct evidence. 
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A. The maximum likelihood technique 

In this section we outline the maximum likelihood procedure used to determine the velocity 
ellipsoid of the Galactic halo. 

Halo stars are assumed to have a gaussian velocity distribution whose principal axes align with 
spherical or cylindrical coordinates. In this appendix we assume the velocity of the Sun relative 
to the galactic center is known and subtracted from stars' radial velocities, and that every star's 
position in space is known exactly; see §2 for a discussion of magnitude and velocity errors. We 
define x as the normalized heliocentric line-of-sight vector to the star, so that the line-of-sight 
velocity vi^g = x • v where v is the true space velocity of the star. The three components of each of 
the vectors x, v, and cr should be thought of as projections along the principal axes of the velocity 
ellipsoid at the specific point in space. The probability (j)dvios of observing a given star with a 
radial velocity between vios and dvios is expressed as an integral over the possible combination of 
true space velocities that would result in an observation between vios and dvios- 

/ (f>l{vi)(l)2{v2)(t)3{v3)dvidv2dV3 (Al) 

-oo J ~oo 

where fi, V2, and ^3 are the components of the star's velocity in the directions of the principal axes. 
Here (j)i{vi) dv = {2'naf)~^l'^ exp {—vf/2af) dv where is the component of the velocity ellipsoid in 
the i direction. Geometry dictates a relation between dfios and the other infinitesimal quanitites: 
dv\os = \xi\dvi. Also, the vector v is constrained by v\os by the relation v\os = x • v so finally our 
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integral becomes 



(/>(vios) dvio 



dv2dV3 4>i 



— {v\o 
Xi 



X2V2 - X3V3) (j)2{v2)4>3{v3 



(A2) 



Here the integral integrates over the 2 and 3 directions, which is mathematically correct if xi 7^ 
0, but to avoid small-number numerical problems, in practice the code shuffles the indices as 
appropriate to ensure that the largest component of x is placed in the denominator. The value of 
the integral is 
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Note that at a given position, stars are most likely to be found at zero velocity (with respect to 
the galactic center) but have a characteristic velocity dispersion cjios, which is a function of the 
geometry of the stars with respect to the Sun and of a. Since the probability of observing a 
set of stars with known space locations and observed line-of-sight velocities is the product of the 
individual probabilities, the log-likelihood is 



N 



4A' 



1 



C']vL-^ln{xlAA') 



(A4) 



where the summation is over all observed stars. Assuming a uniform prior (Lupton 1993), we find 
the maximum of £(o"i, ct2, CJ3) to determine the most likely values of the velocity ellipsoid a. 

If quantities such as the velocity of the Sun with respect to the Galactic center or the bulk 
rotation velocity of the halo are free parameters, they can be determined from this maximum- 
likelihood estimate. These quantities affect vyos for every star, so it is straightforward to extrapolate 
equation A4 to £ = £(<t, v©, fhaio)- 
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Fig. 1. — An illustration of the technique used to incorporate magnitude errors into the maximum 
likelihood analysis, for two artificial stars. The two stars are at the same measured position: 
= (0,30), g = 15.5, and absolute magnitude 0.7. Magnitude errors mean that the stars could 
be located anywhere along the line of sight indicated in the left panel. One star has uios = Okms~^ 
and the other has fios = lOOkms"^. Here we assign a model with <t = (150, 50, 100) km s~^ in 
spherical coordinates, and neglect any velocity of the Sun with respect to the galactic center. The 
upper right panel plots the probability density of the existence of each star with given vy^s at the 
position given by the abscissa (see equation A3). The solid line represents the stationary star, and 
the dotted line represents the receding star. The middle panel shows the weights of Gauss- Hermite 
integration Wj (bottom line) and the Malmquist factor M(r) (top line), as described in the text, 
which are then multiplied to (j) to arrive at the solid (stationary star) and dotted (receding star) 
lines in the bottom panel. The values of (pWjM are summed over the eleven abscissae, and plotted 
in the upper right panel as the horizontal solid and dotted lines, for the stationary and receding 
stars respectively. 
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Fig. 2. — Approximate confidence regions for tfie components of the velocity ellipsoid determined 
from the complete BHB sample (1170 stars) in spherical coordinates, demonstrating Ur is very well 
constrained and that each component is anticorrelated with the other two. In the lower left panel, 
the confidence regions for the halo rotation and the solar velocity in the +?/ direction are shown. 
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Fig. 3. — Each dot represents the resuhs from a Monte Carlo simulation with input cr = 
(150, 100, 80) km s~^ in spherical coordinates. Black dots assume no magnitude or velocity er- 
rors. Blue dots assume the stars are not really at the heliocentric radial positions observed, with an 
error Am = 0.5. Red dots assume the positions are accurate, but that the velocities are measured 
with an error Av = 30kms~^. Green dots assume magnitude and velocity errors, as well as a 10% 
contamination from disk stars. This plot provides evidence that the error estimate as shown in 
figure 2 is accurate. Additionally, the radial component of the velocity ellipsoid is very well con- 
strained, as in the analysis of the real data. One can also make out slight anticorrelations between 
all three parameters, as with the real data. Magnitude errors (blue dots) have no perceivable effect 
on the parameter determinations of the velocity ellipsoid. Velocity errors of 30kms~^ (red dots) 
introduce a small positive shift of order ~ 5kms~^ to each of the three components of the velocity 
ellipsoid, which is most easily seen in the radial component. Contamination (green dots) shifts ag 
and cTr to lesser values and cr,^ to larger values, as expected; also see Table 3. 
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Table 1. Recent results for the local halo velocity ellipsoid. The velocity ellipsoid components 

{R,(j),z) correspond to galactocentric cylindrical coordinates and point radially (outwards), 
retrograde to galactic rotation, and towards the north galactic pole, respectively. Entries in the 
fourth column assume vlsr = (0,220,0) km s-i and V0 — Vlsr = (10, 5, 7) kms ^- 









^^halo 


comments 


source 


125 ± 11 


96 ±9 


88 ±7 




[Fe/H] < -1.2; V <n mag 


Norris, Bessell, & Pickles (1985) 


131 ±6 


106 ±6 


85 ±4 


32± 10 


[Fe/H] < -1.2 


Norris (1986) 


154 ± 18 


102 ± 27 


107 ± 15 


7 ±23 


[Fe/H] < —1.5; red giants 


Carney & Latham (1986) 


133 ±8 


98± 13 


94 ±6 


21 ± 15 


D ~ 4kpc, [Fe/H] < -1.6 


Morrison, Flynn, & Freeman (1990) 


145 ± 10 


100 ± 10 




30± 10 


[Fe/H] < —1.5; kin. sel. subdwarfs 


Ryan & Norris (1991) 






54 ±4 




metal-poor RGs, D < 4kpc 


Ratnatunga & Yoss (1991) 


168 ± 13 


102 ±8 


97 ± 7 


15± 12 


RR Lyr: stat-7r & vios 


Layden et al. (1996) 


151 ±9 


121 ± 7 


103 ±6 


13± 10 


\Z\ < 4kpc, [Fe/H] < -2.2 


Chiba & Beers (2000) 


162 ± 1 


106 ±2 


89 ±2 




Stat vr, {V) -Vq = -216.6 km s'^ 


Gould (2003) 



Table 2. Results 



criteria 


-^stars 


c.s. 


errors 








^halo 


(Jr fsph), (TR (cyl) 


crn (sph), crj, (cyl) 


<jj, (sph), az (cyl) 


all 


1170 


sph 





10.0 


225.0 


7.0 


0.0 


101.4 ± 2.8(2.5) 


97.7 ± 16.4(14.5) 


107.4 ± 16.5(13.9) 


all 


1170 


sph 





10.0 


225.0 


7.0 


-6.0 ± 10.2(10.2) 


101.4 ± 2.8(2.5) 


97.9 ± 16.6(14.5) 


107.9 ± 16.6(13.9) 


all 


1170 


sph 


4.8 


1 ± 6.8(6.5)245.9 ± 13.6(6.7) -7.2 ± 4.0(3.8 


)23.8 ± 20.2(10.1) 


100.6 ± 2.8(2.4) 


94.8 ± 16.6(14.6) 


108.2 ± 16.5(13.8) 


all 


1170 


sph 


m 


10.0 


225.0 


7.0 


0.0 


101.5 ± 2.8(2.5) 


102.1 ± 18.0(15.1) 


103.5 ± 18.2(14.7) 


all 


1170 


sph 


m/v 


10.0 


225.0 


7.0 


0.0 


96.9 ± 3.0(2.6) 


97.5 ± 18.8(15.8) 


99.1 ± 19.0(15.3) 


all 


1170 


cyl 





10.0 


225.0 


7.0 


0.0 


98.7 ± 7.7(5.7) 


104.7 ± 16.8(14.1) 


102.8 ± 5.2(3.3) 


all 


1170 


cyl 


m 


10.0 


225.0 


7.0 


0.0 


99.2 ± 7.8(5.7) 


102.1 ± 17.9(14.8) 


102.9 ± 5.2(3.3) 


all 


1170 


cyl 


m/v 


10.0 


225.0 


7.0 


0.0 


94.5 ± 8.1(6.0) 


97.6 ± 18.7(15.5) 


98.5 ± 5.4(3.5) 


9 < 18 


733 


sph 





10.0 


225.0 


7.0 


0.0 


99.4 ± 4.3(3.4) 


100.0 ± 17.2(14.7) 


110.5 ± 17.1(13.7) 


3 < 18 


733 


cyl 





10.0 


225.0 


7.0 


0.0 


97.2 ± 9.6(7.6) 


109.1 ± 18.0(13.8) 


100.9 ± 6.5(4.3) 


9 < 16 


227 


sph 





10.0 


225.0 


7.0 


0.0 


105.4 ± 10.9(8.2) 


85.7 ± 16.9(14.6) 


117.6 ± 24.5(18.2) 


a < 16 


227 


cyl 





10.0 


225.0 


7.0 


0.0 


108.8 ± 19.7(16.2) 


122.0 ± 23.8(17.5) 


94.1 ± 12.5(8.3) 


9 > 18 


437 


sph 





10.0 


225.0 


7.0 


0.0 


103.2 ± 6.6(3.6) 


28.6 ± 1013.4(741.5) 113.8 ± 211.2(140.1) 


g > 18 


437 


cyl 





10.0 


225.0 


7.0 


0.0 


99.7 ± 13.4(8.6) 


113.3 ± 207.5(142.0) 


104.5 ± 10.4(5.5) 


excl. Sgr 


775 


sph 





10.0 


225.0 


7.0 


0.0 


100.6 ± 3.6(3.0) 


89.9 ± 34.2(26.9) 


106.8 ± 19.5(14.7) 


excl. Sgr 


775 


cyl 





10.0 


225.0 


7.0 


0.0 


96.3 ± 9.5(7.1) 


102.1 ± 18.8(15.2) 


102.6 ± 6.3(4.1) 



Note. — The criteria used to select subsamples from the entire sample of 1170 BHB stars are listed in the first column. The next three columns indicate 
the number of stars in the subsample, the coordinate system used, and whether magnitude or velocity errors were considered in the maximum likelihood 
analysis. The quoted errors not in parentheses are considering correlations among the parameters; the errors in parentheses do not consider correlations. 
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Table 3. Mean and rms of determinations for cr over 100 Monte Carlo simulations for four cases 



error 


color in figure 3 




mean 






rms 






ere 




(Jr 


ere 







black 


150.2 


98.0 


73.7 


3.4 


19.7 


35.6 


m 


blue 


149.5 


98.7 


73.7 


3.9 


18.3 


35.4 


V 


red 


152.3 


101.1 


79.3 


3.8 


19.5 


35.9 


m/v/c 


green 


146.6 


96.8 


104.7 


4.0 


20.5 


30.7 



Note. — The type of error(s) considered is indicated in the leftmost 
column: "0," "m," "v," and "c" refer to no error, magnitude error, velocity 
error, and contamination, respectively. Compare the mean values of cr to 
the input to the simulation of (150, 100, 80) kms~^. 



